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Abstract 

We study theoretically the minimal conductivity of monolayer graphene in the 
presence of Rashba spin-orbit coupling. The Rashba spin-orbit interaction 
causes the low-energy bands to undergo trigonal-warping deformation and for 
energies smaller than the Lifshitz energy, the Fermi circle breaks up into parts, 
forming four separate Dirac cones. We calculate the minimal conductivity for 
an ideal strip of length L and width W within the Landauer-Biittiker formal¬ 
ism in a continuum and in a tight binding model. We show that the minimal 
conductivity depends on the relative orientation of the sample and the probing 
electrodes due to the interference of states related to different Dirac cones. We 
also explore the effects of finite system size and find that the minimal conduc¬ 
tivity can be lowered compared to that of an infinitely wide sample. 

Keywords: mesoscopic systems, quantum wires, carbon 
nanostructnres, charge- and spin- transport 


1. Introdnction 

More than half a century has passed since Landauer derived a formula for 
the conductance of two terminal coherent devices [i|. Then 25 years ago Markus 
Biittiker realized that the two terminal Landauer formula can be extended to 
multi-terminal devices 0. Now, in the literature this approach is commonly 
called Landauer-Biittiker formalism. Over the years it become the standard 
tool for investigating various quantum systems in nanophysics (for a review 
see Refs. saa@). This approach has become an integral part of theoreti¬ 
cal investigations of modern solid states systems such as graphene Q- In the 
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last decade different types of graphene nanostructures proved to be one of the 
most technologically promising and theoretically intriguing solid state systems. 
The dynamics of low energy excitations in graphene is governed by an effective 
Hamiltonian corresponding to massless two dimensional Dirac fermions. Hence 
many physical quantities such as the conductivity, the quantized Hall response 
and optical properties are markedly different from those of conventional two 
dimensional electron systems Q. In bilayer graphene, the interlayer hopping 
results in a trigonally warped Fermi surface which breaks up into four separate 
Dirac cone at low energies. The signatures of this novel electronic structure has 
been studied first experimentally by Novoselov et al. @ and theoretically by 
McCann and Fal’ko (^ . 

Graphene samples, despite the vanishing density of states, show a finite con¬ 
ductivity at the charge neutrality point (at zero Fermi energy). This feature 
of massless Dirac fermions, referred to as minimal conductivity, was intensively 
studied with the Landauer-Biittiker formalism 00, [11. An alternative ap¬ 
proach based on the Kubo formula has also been^plied to study this phe¬ 


nomenon in both monolayer and bilayer graphene IJ, . It was shown that 


in monolayer graphene for wide and short junction the value of the minimal 

4 . ' 


conductivity is (Tq = bilayer graphene neglecting trigonal 

warping the conductivity is tr = 2cro, while including splitting of the Dirac cone 
due to trigonal warping gives extra contributions to the conductivity, increasing 
it to (7 = 6tTo 13 • Later, for finite size of bilayer graphene it was shown by 


Moghaddam and Zareyan m that the trigonal warping results in an anisotropic 
behavior of the minimal conductivity. 

Rashba spin-orbit (RSO) interaction arises once the mirror symmetry of the 
bulk graphene sample is broken by the substrate or an applied electric field per¬ 
pendicular to the graphene sheet. The strength A of the RSO coupling is propor¬ 
tional to this electric field. Photoemission experiments on graphene/Au/Ni(lll) 
heterostructure revealed A ~ 4 meV 0 . Recently, a strong Rashba effect 
with spin-orbit splitting of 70 meA7 has also been observed for graphene on 
Fe(llO) 0. 

Enhanced RSO interaction has a major impact on the transport properties 
of graphene derived samples. Recently the transfer matrix method has been 
employed to study spin dependent transport properties of monolayer graphene 
in the presence of inhomogeneous RSO coupling 0 0. An important conse¬ 
quence of the RSO interaction is that the low-energy behavior of electrons in 
monolayer graphene with RSO coupling is related to that of bilayer graphene 
with trigonal warping but without RSO interaction 2l|. Therefore, we expect 
that the minimal conductivity of monolayer graphene with RSO interaction 
shows a similar anisotropic behavior as that obtained for bilayer graphene in 
Ref. 


161 . 


To see this anisotropic behavior, we calculate the minimal conductivity us¬ 
ing tight binding (TB) calculations and compare it to results obtained from a 
continuous model. We study the effects of finite sample sizes and the crystallo¬ 
graphic orientation as well as the length dependent oscillatory behavior of the 
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Figure 1: Geometry of a graphene device of length L and width W between two electrodes 
doped by potential Uq. Electrons incoming from the left lead are reflected with amplitudes 
r and transmitted with amplitudes t. Between the two contacts we depict the real space 
structure of the monolayer graphene flake (left side) and the energy contours in reciprocal 
space around the K point. The zig-zag direction of the graphene flake makes an angle (/? with 
the electrode interface (y direction). 

minimal conductivity. In our two-terminal calculations, the ballistic scattering 
region of monolayer graphene with length L and width W is contacted by two 
highly doped regions oriented at angle if with respect to the zig-zag direction 
of the graphene lattice (see Fig. [T]). Doping in the electrodes is achieved by 
shifting the Fermi energy with a large potential ffp it is commonly done in 
the literature (see, e.g., Ref. [HI). 

2. Landauer Biittiker formalism for calculating the conductivity 

In the Landauer-Biittiker approach the conductance of a sample is given by 
the transmission probabilities of an electron passing through it: 



( 1 ) 


where tmn are the transmission amplitudes between the propagating modes n 
and m of the left and right electrodes. In what follows, we calculate the minimal 
conductivity in the TB model (for finite W) and compare the results to that 
obtained in the continuous model (for W —>■ oo). Both in TB and continuous 
model the transmission amplitudes tmn are calculated by solving the scattering 
problem of the system. Then the minimal conductivity is defined as cr = -^ G, 
with the conductance G calculated from Eq. m at the charge neutral point of 
graphene, ie, at Ep = 0. 
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Figure 2: Geometry of a graphene sheet. The unit vectors of the hexagonal lattice are ai and 
a 2 , while di = (a 2 — 2ai)/3, d 2 = (ai — 2a2)/3 and da = (ai +a2)/3 are vectors pointing to 
the neighboring atoms. 


2.1. Tight binding model of graphene including RSO coupling 

In the TB model the Hamiltonian Htb of monolayer graphene with RSO 
coupling can be written as 22 |, [21 1 


— Hq 
Ho = -7 

HR = i\ 


Hr, where 

( 2 a) 

E + h.c.) , 

( 2 b) 

T. ^ ~ h.c. . 

( 2 c) 


Here Hq is the Hamiltonian of bulk graphene sheet taking into account only 
nearest neighbor hopping, with hopping amplitude 7 . The operator (uicr) 
creates (annihilates) an electron in the Hh unit cell with spin a on sublattice A, 
while (bja) has the same effect on sublattice B and h.c. stands for hermitian 
conjugate. The unit cell is given by the unit vectors ai and a 2 as shown in 
Fig. [2] The Hamiltonian Hr describes the Rashba spin-orbit interaction where 
s = {sx, Sy, Sz) are the Pauli matrices representing the electron spin, and ijl,v = 
1,2 denote the matrix elements of the Pauli matrices. Here vectors 
connect the nearest neighbor atoms (i, j) pointing from j to i as shown in Fig. [21 


and d is the distance between them, and = d^j ^yd are unit vectors. 

The strength of the spin-orbit coupling is denoted by A which may arise due 
to a perpendicular electric field or interaction with a substrate. 

Using the standard Green’s function techniques | 


24, 


based on the 

Landauer-Biittiker approach we calculate the transmission amplitudes for arm¬ 
chair and zig-zag orientation of the sample. 


2.2. Continuous model of graphene including RSO coupling 

The Hamiltonian of the continuous model as a long wave approximation of 
the TB Hamiltonian Htb in Eq. ([2]) describes low energy excitations around the 
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K and K' points. In our previous publication 12ill w e showed that starting from 
the tight-binding Hamiltonian suggested in Ref. [23 to describe RSO coupling in 
monolayer graphene one can arrive at a form of the Hamiltonian that is unitary 
equivalent to that of bilayer graphene without RSO interaction but including 
the trigonal warping effect due to interlayer hopping l3, 151. In the continuous 
model the Hamiltonian Hk at the K point of the Brillouin zone (BZ) reads as: 


/ 0 

vfP- 

0 

'V\P+\ 
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(3) 


where vp = 3jd/{2h), v\ = 3Xd/{2h), p± = ± ipy and Px,Py are mo¬ 

mentum operators. The Hamiltonian Hp in Eq. © is written in the basis 
(|H t)j \B t)) 1^ i)) \B where {t,4'} refer to spin orientations. A unitary 
equivalent result can be obtained around the Dirac point Kb The four eigen¬ 
values of the Hamiltonian ([^ as a function of the wave number k = [kx, ky) = 
fc(cos a, sin are given by 


E^ik) = ±hvF [kl + (2 + / 12 ) + (-l)n^' 


with 


T = k'l + 2k^ki{2 - -f r/3^(4 -f - Sr A:A^sin(3a), 


(4a) 

(4b) 


where P = v\/vf = A /7 is the dimensionless strength of the spin-orbit coupling, 
k\ = 2/3/c? and n = 1, 2. 

Figure [3] shows the contour plot of the positive and low-energy band and 
the spectrum along the ky direction. The spectrum has a threefold symmetry 
similar to that of bilayer graphene. At moderate energy, direct hopping between 
and leads to trigonal warping of the constant energy lines about each 
valley, but at an energy E less than the Lifshitz energy Ep = 'yP'^ I {A 4 - /3^) the 
effect of trigonal warping is dramatic. It leads to a Lifshitz transition [2^: the 
constant energy line is broken into four pockets, which we refer to as central and 
three leg parts. The Fermi surface is approximately triangle like in the central 
part and each leg part it is elliptical. The distances of the center of the leg parts 
from the K point are kso = 2/3^/c? (see Fig. [^. 

As it has been shown in our previous work [l5j | the Lifshitz transition strongly 
affects the transport properties of monolayer graphene as well. The anisotropy 
of the minimal conductivity in bilayer graphene related to the interference effects 
between the leg parts was predicted by Moghaddam et al. [1^. Therefore, in 
monolayer graphene including the RSO interaction, we also expect a strong 
anisotropy in its conductivity depending on the orientation of the leg parts with 
respect to the electrodes. To see this we calculate the transmission probabilities 
tmn in Eq. ([T]) and the minimal conductivity by solving the scattering problem. 
If we consider the short and wide junction limit {W ^ L), then the electronic 
states can be specified by their energy e and the transverse wavenumber q which 
are conserved during the scattering process. For a given e and q there are four 
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Figure 3: (a) Contour plot of the positive and low-energy band E'^ (in units of A) around 
the K point for {3 = 0.034. Wave vector components kx, ky are in units of "iX/{hvp). The 
center of the pockets has q.I'k jZ rotational symmetry. The distance between the center of the 
pockets and the central Dirac points (points K) is k^o given in the text, (b) The four energy 
bands along the direction ky with kx = 0. 


solutions for the longitudinal wave vector which satisfies the characteristic 
equation det — he] = 0 , where is the 4x4 identity matrix. 

Electronic states in the scattering region {0 < x < L) are denoted by = 

where $ 3 ^, satisfy relation for all possible quantum numbers 1: 

= (5) 


The scattering state between the electrodes is then a linear combination of these 
four electronic states. The longitudinal wave numbers ^be correspond¬ 
ing electronic states in the left (L) and right (R) 

leads can be obtained analogously with a substitution e ^ e — Uq (here Uq is 
the potential on the left and right electrodes as indicated in Fig.[T] If we assume 
an incident states in the L electrode, than the resulted scattering state can be 
written in the form: 


^”(g) 


n’"iq)+ E rn^niqm 

n' — l 

E ^P^SC (9); 

1=1 

E tmniq)'i’R'^iq), 

V. m—1 


(q), X < 0 , 

0 < X < L, 
L < X, 


( 6 ) 


where rn'n{q) and tj^niq) are the reflection and transmission amplitudes, and we 
introduced the arrow —>• (^) to label the right (left) propagating electron states 
in the leads. The reflection and transmission amplitudes have to be determined 
(together with coefficients An) by imposing the continuity condition of the wave 
functions at the interfaces x = 0 and x = L. 
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Figure 4: (Color online) The conductivity (in units of (Tq) of the junction as a function of 
length L (in units of Iso) for (a) zigzag and (b) armchair orientation in continuous model 
obtained from Eq. Q (red solid lines) and from TB calculation (blue dash-dotted and black 
dashed lines) with two different aspect ratios = W/L = 4.71 and = W/L = 6.74 
for zigzag orientation, and R^ = W/L = 3.12 and R^ = W/L = 5.80 for armchair orienta¬ 
tion. The two horizontal dashed lines in Fig. b represent the upper limit of the conductivity 
calculated from the continuous (Sag) and TB model ag) as described in the text. 


Finally, inserting the transmission probabilities tmn('z) into Eq. ([1} we find 
the conductance G. The summation over the transverse wave numbers is re¬ 
placed in a good approximation by the integration ^ f dg. Then the minimal 
conductivity reads: 


c^ = 2 -^G=^L J dqJ2\tmn{q)\^ , (7) 

^ m,n 
— oo ’ 

where in the first equation the factor 2 corresponds to the valley degeneracy. 

3. Results: the minimal conductivity of monolayer graphene with 
RSO interaction 

The minimal conductivity as function of L obtained from the continuous 
model and from TB calculations for zig-zag and armchair orientation are shown 
in Fig. m 

As described in Ref. [l^ the RSO interaction can be characterized by a length 
scale Iso = ti/^SO 1/A^. For short junctions or at low RSO coupling A, that 
is in the limit L/lso 0 , the conductivity for both the armchair and zig-zag 
orientation starts with a{L/lso =0) = ctq. Increasing L/lso the conductivity 
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calculated from the continuous model tends to cr = 3 (To and a = 7j'iaQ for the 
armchair and zig-zag orientation, respectively. 

In the TB calculation for zig-zag orientation, depicted in Fig. 0^, the con¬ 
ductivity closely follows that of the continuous model and tends towards utb = 
7 / 3(70 for longer junctions. Increasing the W/L ratio the subtle peaks of the 
TB and continuous models approach each other. On the other hand, for the 
armchair orientation, shown in Fig.| 4 ) 3 , the results of the TB calculation and the 
continuous model start to deviate for L/lso ~ l-l? that is for increased RSO 
coupling A, tending to a markedly lower value (Ttb = 5/2 (Tq. We also observe 
an enhanced oscillatory behavior as the function of L/lso as compared to the 
calculation done in the zig-zag direction. 

To understand this behavior of the conductivity it is instructive to consider 
the orientation of the Fermi surface around the K and K'points with respect 
to the direction of propagation as shown in Fig. [S] As noted before, trigonal 
warping due to the RSO interaction brakes the Fermi surface into a central 
pocket (dots in Figs. [5Ji and [S^ at the K and K'points) and three extra leg 
pockets labeled by Pi, P 2 and P 3 . 


Zig-Zag edge 

(p=90“ 



Amichair edge 

q )=0 



Figure 5: Schematic drawing of (a) zig-zag and (b) armchair edges. Empty circles and dashed 
lines correspond to the nearest missing sites and bonds to the edge of the ribbon, (c) The 
orientation of one pocket (an ellipse) is given by the angle 0 between the propagation direction 
(k axis) and the semi-major axis of the ellipse, (d) Orientation of the pockets Pi, P 2 , and P 3 
for zig-zag and (e) for armchair edges in the Brillouin zone. 


First we explain the oscillatory behavior of the conductivity shown in Fig. 0)3 
for armchair orientation. In Fig. [5^ the zero energy modes both around the K 
and K' point are at the center of pocket Pi, P 2 and P 3 , and at the center of 
the isotropic Dirac cone. Out of these four modes two (the central Dirac cone 
and pocket P 3 ) have a wave number k = 0 (along the propagating direction) 
and for the other two modes the wave numbers are k = (the centers 

of pocket Pi and P 2 in Fig. 0^. The latter two non-zero propagating modes 
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explain the oscillatory behavior of the conductivity shown in Fig.|3]D. The phase 
shift between the finite k propagating modes (accumulated over one period) for 
an electron bouncing between the electrodes is = ini'/SksoL, where 

TO = 0,1, 2. Then the shortest period of the conductivity is given by A$i = 27r 
from which one finds 


L 

ho 



( 8 ) 


This periodicity can be clearly seen in Fig.|4}D for both TB and continuous cases. 

In the case of the zig-zag orientation all pockets are centered at finite k. In 
both valleys the centers of the central pocket and pocket P 3 are separated by 
kso- this gives L/ho = 2 as the shortest modulation period in agreement with 
our data presented in Fig. 0^. 

Now we explain the marked discrepancy between the continuous model and 
the TB calculations performed in the armchair orientation. For strong RSO 
interaction the conductivity calculated in the tight binding approach can be 
estimated as follows. In general each pocket Pi,P 2 and P 3 shown in Fig. [S] 
corresponds to one anisotropic Dirac cone and gives a separate contribution to 
the conductivity. The total conductivity is given by 


cr = ncCTc + ^where 


r(0.) = 


vl cos^( 0 i) -b sin^( 0 i) ctq 
VaVb 4 ’ 


(9a) 

(9b) 


and ac = c^q/A is the contribution from the central Dirac cone ll| while (j(0i) 
is the minimal conductivity related to a single anisotropic Dirac cone. This 
result was first derived by Nilsson et al. in Ref. 27|. Note that the same result 
can be obtained by the general approach developed in Ref. [l^ . Here nc and rii 
are the number of open channels for the central Dirac cone and the leg pocket 
Pi, respectively, and Va and Vf, are the Fermi velocities along the two principal 
axes of the ellipse corresponding to the pocket Pi with i = 1,2,3. For our case 
Vb = 3va [27| for the three legs. 0^ is the angle of the direction of the semi-major 


axis of the ellipse with respect to the direction of propagation (see Fig. [S]:). 
One can see from Fig. [5] that around the K point for armchair orientation 
0i = Ttt/G; IItt/G; 7 r /2 for pocket Pi, P2 and P3, respectively, while for zig¬ 
zag orientation 0i = Stt/S; tt/S; tt around the K point, and 0i = 47 r/ 3 ; 27r/3; 0 
around the K' point for the pocket Pi, P2 and P3, respectively (see Fig. B- 
We now determine the number of open channels nc and Ui in Eq. (jH]). The 
boundary condition for the system demands that the wave function at the two 
edges of the ribbons should be zero at the empty sites shown schematically in 
Figs. [S^ andjSJi. Thus, including the spin we have four equations to satisfy the 
boundary conditions. For a given energy e and wave number k corresponding 
to the propagating mode the possible transverse modes can be calculated from 
the dispersion relation. Graphically it means that these transverse modes can 
be obtained by drawing a vertical line at a given k that intersects the given 
constant energy e contour. For example, for zig-zag orientation for wave number 
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k for which the vertical line passes through the center of pocket P3 around the K 
point there are two transverse modes (the vertical line crosses the energy contour 
at two points in Fig. [^i, while for wave number k for which the vertical line 
passes through the center of pockets Pi and P2 we have four transverse modes. 
Hence, it follows that in the first case the number of open channels ni = 0 
since the four boundary conditions cannot be satisfied by two transverse modes. 
Similarly, for the central isotropic Dirac cone ric = 0. However, for the second 
case ni = 712 = 1 because we have four transverse modes. The same is true for 
the propagating mode k around the K' point (valley degeneracy). In summary, 
the open channels for zig-zag ribbons are nc = 0 and ni = n2 = 2,113 = 0- 
From a similar consideration we find that for armchair orientation nc = 2 and 
111 = 712 = = 2. Thus the minimal conductivity of monolayer graphene 

with RSO coupling given by ([9]) is cr = 7/3 cto for zig-zag and a = for the 

armchair orientation, in very good agreement with the TB calculations. 

4. Conclusions 

We have investigated the minimal conductivity of monolayer graphene in 
the presence of Rashba spin-orbit interaction. We have employed tight binding 
calculations and a continuous model, to determine the interplay of the crystal¬ 
lographic orientation of the sample with the anisotropic nature of the minimal 
conductivity. Contrasting the results obtained for a graphene strip of finite 
width to that of an infinitely wide sample, we show that the boundary condi¬ 
tion for a finite flake may, depending on the orientation, reduce the value of 
the minimal conductivity compared to that of the infinitely wide. All our cal¬ 
culations have been performed in the spirit of the Landauer-Biittiker approach. 
We hope our results are a tribute for the long lasting legacy of this simple yet 
powerful formalism and to the memory of Markus Biittiker. 
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